Methods for Solving the Traveling Salesman 



Technical Field 



The present invention relates to computational models, simulation, and other 
techniques for solving complex problems in combinatorics and artificial intelligence and, 
in particular, it relates to methods for solving the Traveling Salesman Problem, whxch 
have a broad range of applications. 



Background 



In the traditional Traveling Salesman Problem (TSP), a traveling salesman must 
visit customers in a number of cities and he knows the distance between pairs of cues he 
wants to visit. Since he is interested in not being too long on the road, he wants to taie 
the shortest tour. FIG. 1 shows a number of cities to be visited in a tradthonal TSP. 
Suppose a traveling salesman must visit the following 8 capitol cities shown in FIG. 1: 
(a) Boise, ID; (b) Salt Lake City, UT; (c) Phoenix, AZ; <d) St. Paul, MN; (e) Columbus, 
OH- (0 Jefferson City, MO; (g) Richmond, VA; and (h) Austin, TX. FIG. 2 shows the 
shortest tour between the cities shown in FIG. 1 starting with a (Boise, ID) and endtng 
with a (Boise, ID). The shortest tour is represented by the ordering a. b, c. h, f, g. e, d, 
which means he starts in city a (Boise, ID), then goes to cities b. c, h, f, g, e, and d m that 
order (Salt Lake City, UT to Phoenix, AZ to Austin, TX to Jefferson City, MO to 
Richmond, VA to Columbus, OH to St. Paul, MN) and returning to city a (Bo.se, ID). 
However, for various real-world applications of the TSP, there is a need for a modtfied 
solution to the TSP mat starts and ends in different cities. FIG. 3 shows an alternate 
solution of the TSP shown in FIG. 1 in which the tour a, b, c, h,f, d, e, g begins and ends 
in different cities. 

Furthermore, the tour shown in FIG. 2 is an idealized solution. FIG. 4 shows a 
more realistic example solution of the TSP than the example solution shown in FIG. 2. 
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To travel by road would require a more roundabout path. For example, one would 
probably take highway 84 East from Boise, ID to Sal. Lake City, highways 15 South, 93 
South and 40 East to Phoenix, AZ, highways 17 North, 40 East, 287 Sou*, and 35 Sou* 
to Austin, TX, back on highway 35 North to highways 75 North 69 North , 44 East , 73 
North and 54 North to Jefferson City, MO, highways 54 North, 70 East, and 64 East to 
Richmond, V A, back on highway 64 West to highways 8 1 North and 70 West to 
Columbus, OH, highways 70 West , 65 North, 80 West, 39 North, 90 North, and 94 West 
,o S. Paul, MN and highways 94 West, 15 South, 86 West, and 84 West hack to Botse, 
m Also, extra time would be needed on slower roads, such as roads with construct™ or 
roads through large cities a. rush hour. Travel by plane, boa,, or train also has certatn 
routes, congestion, and delays. lh» is a need for a solution to the TSP that takes into 
account the actual or realistic traveling distance, costs and time. 

The TSP is an optimization problem, which is concerned with choosmg the best 
among all possible arrangements, according to some criteria, such as distance, time, and 
traveling expenses. Therefore, it is not sufficient to find an arbitrary solution; instead one 
is interested in the best -or a. least a very good-solution. However, it is no, practical 
to enumerate all possible ordered paths for a large number of cities in order ,o dew*, 
the absolute best one. For a tour of . cities, the enumeration would have «! (n-factonal) 
ordered paths, where »! = » x („ - 1) x („ - 2) x x 1 . The function n-factorial increases 
exponentially. For example, the enumeration of all possible paths for 8 cities ts 8! - 8 x 
7x6x5x4x3x2x1 = 40,320. For 100 cities, the enumeration of all possible pattis is 
,00' which is approximately 9.3 x 10'" ( which is "93" followed by 156 zeros). If a 
memod enumerated one hundred million pattis per second, i.e., a, 100 MIPS, die speed of 
a typical PC, it would take approximately 2.9 x 10'" years to enumerate all of them! (2.9 
x 10'" years is approximately = 9.3 x 10'" paths * (1,000,000 paths/sec x 60 sec/mm x 
60 min/hr x 24 hr/day x 365 days/yr).) Clearly, there is a need to produce a result while 
we are still alive. The TSP has been called an NP-hard problem because no method 
exisK for computing Ore solution in polynomial time. Polynomial time means a computer 
naming time or an order of complexity of OOH or better, where * is a constant. There is 
, a need for a method for solving the TSP in polynomial time. 
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Various methods have tried to solve the TSP. These methods can be broadly 
classified as heuristic methods, natural optimization techniques (such as linear 
programming, mixed integer programming, and quadratic programming), neural 
networks, and techniques based on Boltzman machines, such as simulated annealing. 
Heuristic and natural optimization techniques only give local optima rather than global 
optima, their order of complexity is 0(«!) and they are unstable when a large number of 
cities are considered. Neural networks are dependent on a training methodology. 
Simulated annealing gives a solution, but in doing so it takes a lot of time or is slow at 
arriving at the solution for a large number of cities. There is a need for a method that 
gives a solution in less time for any number of cities, with complexity of 0(n + *), where 
A: is a constant. 

For these and other reasons, there is a need for the present invention. 



Summary 

The present invention is directed to several aspects of methods for solving the 
TSP including the use of a genetic algorithm. Different aspects of the present invention 
include a solution to the TSP that starts and ends in different cities, takes into account 
realistic traveling costs, and completes in polynomial time, such as 0(n + *), where * is a 

constant. 

One aspect of the present invention is a computer readable medium having 
computer-executable instructions for performing a method for solving a traveling 
salesman problem. According to the method, a set of locations to visit are selected. A 
starting point and an ending point from the set of locations is selected. Then, a search 
method using a genetic algorithm is applied to the set of locations and a route is provided 
as a solution to the traveling salesman problem. The genetic algorithm includes an 
objective function that simultaneously minimizes distance and time. The route is from 
the starting point to the ending point, visiting all locations in the set of locations. The 
objective function evaluated the route to a value that is lower than any other path 
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searched. By simultaneously minimizing distance and time, the genetic algorithm takes 
into account realistic traveling costs and is better than evaluating distance alone. 

Another aspect of the present invention is a computer readable medium having 
computer-executable instructions for performing a method for finding a route in a supply 
chain. Finding a route in a supply chain is an example of a solution to the TSP. 
According to the method, a set of orders are accessed. A set of supplier locations are 
selected to visit to pick up products to fill the set of orders. The method searches for 
pickup paths visiting each supplier location in the set of supplier locations. The 
searching is done by applying a genetic algorithm that simultaneously minimizes distance 
and time. The method provides a pickup route as a solution to the TSP. The pickup route 
comprises a pickup path that best simultaneously minimizes distance and time compared 
with all the other pickup paths searched. The genetic algorithm has an order of 
complexity of 0(« + k), where k is a constant, so it completes in polynomial time. 

Another aspect of the present invention is a method of searching among locations 
for paths. Searching among locations for paths is another example of a solution to the 
TSP. According to the method, an initial population of chromosomes is randomly 
generated to be a current population. A fitness for each chromosome in the current 
population is calculated. One or more pairs of parent chromosomes to generate a new 
population are selected. The pairs are crossed over to generate the new population. The 
offspring are mutated to modify the new population. The current population is made a 
previous population and the current population is replaced with the new population. New 
generations are formed by repeating the calculating, selecting, crossing over, mutating, 
and making acts, until a previous best fitness in the previous population is the same as a 
current best fitness in the current population. This method provides a path represented by 
a chromosome having a best fitness in the current population. Each chromosome holds a 
predetermined number of integers. One integer is located at each position in the 
chromosome. Each integer represents one location from a set of locations. The fitness 
represents distance and time simultaneously. The pairs of parent chromosomes are chosen 
from the current population based on their fitness. The pairs are crossed over at a 
randomly chosen point to form offspring. The crossover is carried out with a crossover 
probability p c between about 0.6 and 0.9. Mutation is done at each position on the 
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chromosomes of the offspring. The mutation is carried out with a mutation probability 

p m between about 0.001 and 0.01. 

Another aspect of the present invention is a computer system for finding 
distribution chains to satisfy customer demand. Finding distribution chains to satisfy 
customer demand is another solution to the TSP. The computer system comprises a 
storage device, an output device, and a processor programmed for repeatedly performing 
a method. According to the method, current orders are determined. A set of locations to 
visit which comprises locations of suppliers with products available for the current orders 
is determined. Also, a starting location in the set of locations and a distribution center 
proximate to customer locations are determined. A genetic algorithm is applied to find a 
first route while minimizing distance and time. The method provides the first route on 
the output device, stores any undelivered products in a buffer, and removes the orders 
that the first route will fill from the current orders. The customer locations are delivery 
locations for each customer who placed one of the current orders. The first route starts at 
the starting location, ends at the distribution center, and visits each location in the set of 
locations. The undelivered products are stored in the buffer on the storage device and 
will be included in the current orders during the next iteration of the method. 

Another aspect of the present invention is a data signal having a data structure 
stored thereon, which comprises a first field containing locations and a second field 
containing a request for a route through the locations in the first field. Another aspect of 

the present invention is a data signal having a data structure stored thereon, which 

comprises a field containing a route which comprises ordered locations to be visited. 
Another aspect of the present invention is a method of communication, 

comprising: receiving a plurality of locations and a request for a route through the 

locations, and sending, in response to the request, the route. 

These and other features, aspects, and advantages of the present invention will 

become better understood with reference to the following description and appended 

claims. 
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Brief Descrip *'"" «f the Drawings 

FIG. 1 shows a number of cities to be visited in a traditional Traveling Salesman 
Problem (TSP). 

1 i„*;«r. rt f the- T<?P shown in FIG. 1 in which a tour 
FIG. 2 shows an example solution ot the im* snown 

begins and ends in the same city. 
FIG. 3 shows an alternate solution of the TSP shown in FIG. 1 in which a tour 

begins and ends in different cities. 
FIG. 4 shows a more realistic example solution of the TSP than the example 

solution shown in FIG. 2. 
FIG 5 shows an example solution of the TSP shown in FIG. 1 and an example 
encoding of the solution in a chromosome, according to embodiments of 
the present invention. 
FIG. 6A shows an overview of an example embodiment of a method for solving 
the TSP applied to supply chain management, according to embodiments 
? Z of the present invention. 

S FIG. 6B shows a block diagram of a computerized system to implement methods 

h for solving the TSP, as further described in FIGS. 7-18. 

£ FIG. 7 shows a flow chart of an example embodiment of a method for solving the 

I 2Q TSP to be implemented by the general processing modules of FIG. 6B. 

S FIG. 8 shows a more detailed flow chart of alternate embodiments of the method 

for solving the TSP of FIG. 7. 
FIG. 9 shows a flow chart of an example embodiment of a method for finding a 
route in a supply chain to be implemented by the general processing 
modules of FIG. 6B. 

FIG. 10 shows a more detailed flow chart of alternate embodiments of the method 

for finding a route in a supply chain of FIG. 9. 
FIG. 11 shows another more detailed flow chart of alternate embodiments of the 
method for finding a route in a supply chain of FIG. 9. 
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FIG. 12 shows a flow chart of an example embodiment of a method of searching 
among locations for paths to be implemented by the general processing 

modules of FIG. 6B. 
FIG. 13 shows a more detailed flow chart of alternate embodiments of the method 

of searching among locations for paths of FIG. 12. 
FIG. 14 shows another more detailed flow chart of alternate embodiments of the 

method of searching among locations for paths of FIG. 12. 
FIG. 15 shows a flow chart of an example embodiment of a method for finding 
distribution chains to satisfy customer demand to be programmed on a 
processor in a computer system. 
FIG. 16 shows a more detailed flow chart of alternate embodiments of the method 

for finding distribution chains to satisfy customer demand of FIG. 15. 
FIG. 17 shows example embodiments of data signals to be communicated by an 

example embodiment of a method for solving the TSP. 
FIG. 18 shows an example embodiment of a method of communication for a user 

interface to a method for solving the TSP. 
FIG. 19 shows an example of a suitable computing system environment for 

implementing embodiments of the present invention, such as thbse shown 
in FIGS. 6A,6B, and 7-18. 

Detailed Description 

In the following description, reference is made to the accompanying drawings 
which form a part hereof, and in which is shown by way of illustration specific 
embodiments in which the invention may be practiced. These embodiments are 
described in sufficient detail to enable those skilled in the art to practice the invention, 
and it is to be understood that other embodiments may be utilized and that structural, 
logical and electrical changes may be made without departing from the scope of the 
present invention. The following description is, therefore, not to be taken in a limited 
sense, and the scope of the present invention is defined by the appended claims. 
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The present invention is directed to various aspects of methods for solving the 
Traveling Salesman Problem (TSP). Novel methods for solving the TSP and example 
applications of the methods are described below. The following detailed description 
includes representing location information for a TSP as chromosomes. By representing 
or encoding location information as chromosomes, e.g. storing it in a chromosome data 
structure, methods may incorporate genetic algorithms. The methods for solving the TSP 
may be applied to a broad range of applications, such as supply chain management. 

FIG. 5 shows an example solution 500 of the TSP shown in FIG. 1 and an 
example encoding of the solution in a chromosome, according to embodiments of the 
present invention. In FIG. 5, the eight cities of FIG. 1 are represented by the integers 1-8 
as follows: (1) Boise, ID; (2) Salt Lake City, UT; (3) Phoenix, AZ; (4) St. Paul, MN; (5) 
Columbus, OH; (6) Jefferson City, MO; (7) Richmond, VA; and (8) Austin, TX. In other 
words, the letters, a, b, c, d, e,f,g,h in FIGS. 2-4 have been replaced by the integers 1, 2, 
3, 4, 5, 6, 7 ,8 in FIG. 5. The example solution 500 of the TSP, shown in FIG. 5, may 
also be written as the ordered path 1, 2, 3, 8, 6, 4, 5, 7, which represents a solution 500 of 
the TSP from (1) Boise, ID to (2) Salt Lake City, UT to (3) Phoenix, AX to (8) Austin, 
TX to (6) Jefferson City, MO to (4) St. Paul, MN to (5) Columbus, OH to (7) Richmond, 
VA. 

FIG. 5 also shows how the tour 500 may be represented in a chromosome data 
structure 502. A chromosome data structure 502 in a genetic algorithm refers to a 
candidate solution to the TSP. Some other candidate solutions to the TSP for 8 cities are 
{I, 4, 2, 3, 8, 6, 5, 7}, {2, 1, 6, 8, 3, 4, 5, 7}, {7, 8, 6, 5, 4, 3, 2, 1} and other permutations 
of the numbers 1-8. Once location information is encoded into a chromosome data 
structure, then methods including a genetic algorithm may be used to solve the TSP. 

A genetic algorithm is an evolutionary computation method, which is an 
abstraction or model of biological evolution. In a genetic algorithm, candidate solutions 
are reproduced and evaluated for fitness in successive generations until a best solution is 
found using a kind of "natural selection" together with genetics-inspired operations like 
crossover or replication and mutation. In biology, each gene is located at a particular 
position on the chromosome. In a genetic algorithm solving the TSP, each city may be 
represented as an integer at a particular position in the chromosome data structure. For 
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example, city (1) Boise, ID is encoded in the first position 504 of the chromosome data 
structure 502. Each element of the solution 500, which is the ordered path 1,2,3, 8, 6, 4, 
5, 7, is encoded into the chromosome data structure 502 in its corresponding position, i.e. 
the first element 1 is in the first position, the fourth element 5 is in the fourth position, 
etc. 

Novel methods of solving the TSP will be described in detail below. These novel 
methods can be applied to a broad range of applications. One application that will benefit 
from the new methods for solving the TSP is supply chain management. Supply chain 
management includes timely distribution of goods to meet customer demand. A supply 
chain is a business process that enables customer demand for a product or service to be 
satisfied. Supply chain management decisions include logistics, production, inventory, 
transportation, and scheduling. A method to find solutions to the TSP aids or automates 
supply chain management decisions by, for example, generating cost effective 
distribution paths. One example of supply chain management is a web-based trading 
exchange or portal connecting buyers and sellers along the supply chain of the 
automotive industry. Customer demands for any automotive products in any amount 
from any part of the world need to be satisfied in time. 

FIG. 6A shows an overview of an example embodiment of a method for solving 
the TSP applied to supply chain management, according to embodiments of the present 
invention. One example of a supply chain management application is a typical "Dot Com 
Shop" operation 600 like the one shown in FIG. 6A. A "Dot Com Shop" 602 provides 
online ordering for customers 604 over a network 606. The network 606 may be a local 
area network (LAN), a wide area network (WAN), such as the Internet, an intranet, or 
any other kind of network, for example, a telephone network. A method of solving the 
TSP in the form of processing modules 608 running on a computing device receives the 
orders placed online by customers 604 from the "Dot Com Shop" 602, over the network 
606, or from a database 610. Processing modules 608 are computer-executable 
instructions that perform methods and may be implemented in any software language or 
script or in hardware or firmware in the computing device. 

However, there are many other applications of methods for solving the TSP 
besides supply chain management, including optimization tasks (e.g., numerical 
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optimization, circuit layout, and job scheduling), automatic programming (e.g., sorting 
networks), machine learning (e.g., neural networks, symbolic production systems, and 
sensors for robots), economics (e.g., bidding strategies), immune system models, 
ecological models, population genetics, evolution and learning, and social systems (e.g., 
the evolution of cooperation and communication in multi-agent systems). 

In FIG. 6A, the database 610 may be associated with or in communication with 
the computing device and may contain various kinds of information, such as order 
information, operational constraint information, customer information, supplier 
information, and other information. Order information comprises order identification 
numbers, customer name and address, delivery requirements, name of products, type of 
products, number of products, price of products, total price of order, and the like. 
Operational constraint information comprises customer preferences, supplier or 
warehouse policies and procedures, travel conditions, transportation requirements, 
equipment malfunctions, and the like. Customer information comprises name, address, 
delivery location, delivery requirements, order information and history, and the like. 
Supplier information is information concerning suppliers, such as manufacturing 
operations and warehouses. Supplier information comprises name, address, product 
availability, lines of products, and the like. 

In FIG. 6A, the processing modules may produce solutions to the TSP based on 
input from the database in the form of pickup routes 612 and delivery routes 614. The 
routes may be for any kind of transportation, such as plane, ship, truck, or the like. A 
pickup route 612 is an ordered path, such as St, S 2 , S 3 , S 8 , S 6 , S 4 , S 5 , S 7 , representing a 
plan to pickup products starting with Supplier in Boise, ID, then to Supplier in Salt 
Lake City, UT, to Supplier in Phoenix, AZ, to Suppliers in Austin, TX, to Supplier in 
Jefferson City, MO, to Supplier in St. Paul, MN, to Suppliers in Columbus, OH, and 
ending at Supplier in Richmond, VA. The routes are not limited to the United States, 
but may be anywhere in the World. A delivery route 614 is an ordered path, such as C,, 
C 2 , C 3 , C 4 , representing a plan to deliver products to customers starting with Customer to 
Customer to Customer and ending at Customer^ The pickup route 612 and the delivery 
route 614 may be linked by a distribution center 616. The distribution center 616 may be 
at one of the suppliers or a separate location. The distribution center 616 is a place where 
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products are shipped for local or regional delivery to a set of customers. For example, a 
distribution center may be located in a major metropolitan area and service customers 
throughout the metro area and the surrounding suburbs. Any shipping method may be 
used for delivery, such as vans, trucks, and the like. 

FIG. 6B shows a block diagram of a computerized system to implement methods 
for solving the TSP, as further described in FIGS. 7-18. As shown in FIG. 6B, the 
computerized system includes processing modules 620, inputs 622, and outputs 624. The 
general processing modules 620 in FIG. 6B correspond to the processing modules 608 
running on a computing device in FIG. 6A. Various aspects and embodiments of the 
inputs 622 and outputs 624 will become clear with reference to the following description. 

FIG. 7 shows a flow chart of an example embodiment of a method 700 for solving 
the TSP to be implemented by the general processing modules of FIG. 6B. The method 
700, shown in FIG. 7, may be implemented as processing modules 620, inputs 622, and 
outputs 624, as shown in FIG. 6B. The inputs 622 comprise a set of locations and the 
outputs 624 comprise a route as a solution to the TSP. Depending on the application, 
locations may be cities, locations of shippers, locations of warehouses, locations of 
manufacturers, locations on a circuit board, positions of a robot arm, and other kinds of 
locations. Similarly, a route may be an ordered list of cities, a design for a circuit board, 
a movement for a robot arm, and other kinds of routes. 

As shown in FIG. 7, one aspect of the present invention is a computer readable 
medium having computer-executable instructions for performing a method 700 for 
solving a TSP. The method 700 comprises: selecting a set of locations to visit 702, 
selecting a starting point and an ending point from the set of locations 704, applying a 
search method to the set of locations 706, and providing a route as a solution to the TSP 
708. The route is from the starting point to the ending point, visiting all locations in the 
set of locations. 

An aspect of some embodiments of the present invention is that the search method 
includes a genetic algorithm. In one embodiment, the search method is a combinatoric 
approach to the genetic algorithm. One example embodiment of a search method 
including a genetic algorithm is shown in FIG. 12. The present invention, however, 
comprehends various other embodiments. Combinatorics includes a set of techniques for 
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choosing, among all possible arrangements, that which is best according to some criteria, 
such as matching and covering. 

An example genetic algorithm is shown in Table 1 . Other example embodiments 
of a genetic algorithm may also be used in the present invention. The genetic algorithm 
5 includes an objective function to simultaneously minimize distance and time. The 

objective function would correspond to f(x) in Table 1. An objective function is used to 
evaluate the fitness of each path in the population according to some criteria, such as 
distance and time. A population is a set of paths or chromosomes generated by the 
genetic algorithm. The objective function evaluated the route, which was chosen as the 
10 solution to the TSP, to a value that is lower than any other path searched. In other words, 
the genetic algorithm evaluated each path in each generation and, of all the paths 
evaluated, the one with the best fitness was chosen as the solution to the TSP. This is 
*y because generally the fittest paths or chromosomes survive to replicate the next 

jj= generation. 

W i5 A genetic algorithm typically contains three operators: selection, crossover, and 

*B mutation, as illustrated in Table 1 . Selection is an operator that selects chromosomes in 

1^ the population for reproduction. The fitter the chromosome, the more times it is likely to 

O be selected to reproduce. Crossover is an operator that randomly chooses a position on a 

rb chromosome and exchanges subsequences before and after the position between two 

:L' 20 chromosomes to create offspring. For example, the ordered paths 3, 2, 7, 4, 5, 6, 7, 8 and 
^ 7, 2, 3, 8, 6, 4, 5, 7 could be crossed over after the third position to produce the offspring 

J, 2, 7, 8, 6, 4, 5, 7 and 7, 2, 5, 4, 5, 6, 7, 8. The crossover operator roughly mimics 
biological recombination between two single-chromosome organisms. Crossover is done 
with a probability such as 0.6. Mutation is an operator that randomly flips some of the 
25 positions in a chromosome. For example, the path 3, 2, 1, 4, 5, 6, 7, 8 might be mutated 
in its second and seventh positions to yield 3, 7, 7, 4, 5, 6, 2, 8. Mutation can occur at 
each position in a chromosome with some small probability, such as 0.001. 

Table 1. An Example Genetic Algorithm 

1 . Start with a randomly generated population of n size chromosomes. 

2. Calcu late the fitness f(x) of each chromosome x in the population. 
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3. Repeat the following steps until n offspring have been created: 

a. Select one or more pairs of parent chromosomes from the current population, the 
probability of selection being an increasing function of fitness. Selection is done "with 
replacement" meaning that the same chromosome can be selected more than once to 
become a parent. 

b. With crossover probability p c , cross over the pairs at a randomly chosen point to 
form two offspring. If no crossover takes place, form two offspring that are exact copies 
of their respective parents. 

c. Mutate the two offspring at each position with mutation probability p m and place the 
resulting chromosomes in the new population. If n is odd, one new population member 
can be discarded at random. 

4. Replace the current population with the new population. 

5. Go to step 2. 

FIG. 8 shows a more detailed flow chart of alternate embodiments 800 of the 
method 700 for solving the TSP of FIG. 7. In one embodiment of the method 700, the 
genetic algorithm generates populations of size 2*, where the set of locations is of size k 

5 802. In another embodiment, the genetic algorithm comprises a selection operation, a 
crossover operation, and a mutation operation 804. In another embodiment, the selection 
operation is CHC selection 806. CHC selection is characterized by elitist selection, 
incest prevention, a highly disruptive form of crossover, and restarts. Elitist selection is 
like survival of the fittest. Elitist selection involves cross-generational competition by 

10 merging parent populations and child populations and then ranking then according to 
fitness. The parent populations are created by selecting a predetermined number of the 
fittest members of a population. Incest prevention is a way to slow convergence to a 
solution by not mating parents that are too similar, such as parents with a difference 
threshold of a Hamming distance of l/4 th the number of locations to visit in the route and 

15 by decreasing the difference threshold if no children can be created. Restarts prevent 

premature convergence by reintroducing diversity when a population is nearly converged. 
This is done by performing a mutation operation. In another embodiment, the crossover 
operation is greedy crossover 808. Greedy crossover is characterized by always choosing 
a location in the parent chromosome for the next position in an offspring chromosome 

20 that is closest to the location in the previous position in the offspring chromosome, if 
such a location exists. It is greedy in the sense of attempting to minimize distance. In 
another embodiment, the crossover operation has a crossover probability between about 
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0.6 and 0.9 810. In another embodiment, the mutation operation comprises swapping two 
locations at two random positions in a chromosome 812. In another embodiment, the 
mutation operation has a mutation probability between about 0.001 and 0.01 814. In 
another embodiment, the genetic algorithm further comprises: encoding locations^into 
chromosomes so that the chromosomes have the same number of positions as k and the 
locations are represented in each position by an integer 816. Locations may also be 
represented by bits (e.g. 010001 1 1) or characters, such as a, b, c, d, e,f,g,h in the 
present invention. 

Novel methods for solving the TSP have been described by reference to FIGS. 7 
and 8. As described earlier, solutions to the TSP may be applied to supply chain 
management. FIG. 9 shows a flow chart of an example embodiment of a method for 
finding a route in a supply chain to be implemented by the general processing modules of 
FIG. 6B. The method 900, shown in FIG. 9, may be implemented as processing modules 
620, inputs 622, and outputs 624, as shown in FIG. 6B. The inputs 622 comprise a set of 
orders and a set of supplier locations to visit and the outputs 624 comprise a pickup route. 
One aspect of the present invention is a computer readable medium having computer- 
executable instructions for performing a method for finding a route in a supply chain 900. 
The method 900 comprises: accessing a set of orders 902, selecting a set of supplier 
locations to visit to pick up products to fill the set of orders 904, searching for pickup 
paths visiting each supplier location in the set of supplier locations, by applying a genetic 
algorithm that simultaneously minimizes distance and time 906, and providing a pickup 
route which comprises a pickup path that best simultaneously minimizes distance and 
time compared with all the other pickup paths searched 908. 

FIG. 10 shows a more detailed flow chart of alternate embodiments 1000 of the 
method 900 for finding a route in a supply chain of FIG. 9. In one embodiment, 
accessing a set of orders 902 further comprises accessing current orders and buffered 
orders 1002. Orders that cannot be delivered in the current route are stored in a buffer so 
that they can be delivered in a later route. For example, orders may not be able to be 
filled due to products not being available by any supplier until a later date. Also, orders 
may be continuously be placed by customers while the method 900 is run only 
periodically. In that case, orders placed between runs will be buffered. In another 
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embodiment, accessing a set of orders 902 further comprises accessing orders for a 
particular delivery time period 1004. In another embodiment, the method 900 further 
comprises accessing operational constraint information 1006. In another embodiment, 
the method 900 further comprises satisfying the operational constraint information 1008. 

FIG. 1 1 shows another more detailed flow chart of alternate embodiments of the 
method for finding a route in a supply chain of FIG. 9. In one embodiment, selecting a 
set of supplier locations to visit to pick up products to fill the set of orders 904 further 
comprises computing a distance between each supplier location in the set of supplier 
locations 1 102. In another embodiment, the set of supplier locations includes a 
distribution center 1 104. In another embodiment, selecting a set of supplier locations to 
visit to pick up products to fill the set of orders 904 further comprises eliminating 
redundancies in the set of supplier locations by selecting a first supplier over a second 
supplier, when the first supplier's available products and second supplier's available 
products are the same, but the first supplier has more available products to fill the orders 
than the second supplier 1 106. For example, the method 900 may check whether a 
particular product is available at only one or many suppliers. If it is available at only one 
supplier, then that supplier is selected. Otherwise, the supplier that can supply the largest 
number of those products is selected. In another embodiment, searching for pickup paths 
906 further comprises searching for pickup paths ending at the distribution center 1 108. 
In another embodiment, the pickup route provided in 908 is a pickup route that ends at 
the distribution center 1110. In another embodiment, the method 900 further comprises 
determining a set of customer locations for each customer placing an order in the set of 
orders. In another embodiment, the set of customer locations includes a distribution 
center 1114. In another embodiment, the method 900 further comprises computing a 
distance between each customer location in the set of customer locations 1116. In 
another embodiment, the method 900 further comprises searching for delivery paths 
visiting each customer location in the set of customer locations, by applying a genetic 
algorithm that simultaneously minimizes distance and time 1118. In another 
embodiment, searching for delivery paths includes searching for delivery paths starting at 
the distribution center 1 120. In another embodiment, the method 900 further comprises 
providing a delivery route, which comprises a path that best simultaneously minimizes 
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distance and time compared with all the other delivery paths searched 1 122. In another 
embodiment, the delivery route begins at the distribution center 1 124. 

FIG. 12 shows a flow chart of an example embodiment of a method of searching 
among locations for paths to be implemented by the general processing modules of FIG. 
6B. One aspect of the present invention is a method of searching among locations for 
paths 1200. According to some embodiments of the present invention, the locations are 
encoded as a chromosome data structure. The method 1200 comprises: randomly 
generating an initial population of chromosomes to be a current population, each 
chromosome holding a predetermined number of integers with one integer at each 
position in the chromosome, each integer representing one location from a set of 
locations 1202, calculating a fitness for each chromosome in the current population, the 
fitness representing distance and time simultaneously 1204, selecting one or more pairs 
of parent chromosomes from the current population based on their fitness to generate a 
new population 1206, crossing over the pairs at a randomly chosen point, with crossover 
probability p c between about 0.6 and 0.9, to form offspring to generate the new 
population 1208, mutating the offspring at each position on their chromosomes, with 
mutation probability p m between about 0.001 and 0.01, to modify the new population 
1210, making the current population a previous population and replacing the current 
population with the new population 1212, forming new generations by repeating the 
calculating, selecting, crossing over, mutating, and making acts, until a previous best 
fitness in the previous population is the same as a current best fitness in the current 
population 1214 and providing a path represented by a chromosome having a best fitness 
in the current population 1216. 

An example of a randomly generated initial population of chromosomes 
representing ordered paths through locations 1-8 is shown in Table 2. 

Table 2. An Example of A Randomly Generated Initial Population 

35782146 
43827 561 
68 34 5721 
74 236815 
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An example of calculated fitness values for each chromosome in the population of Table 
2 is shown in Table 3. 

Table 3. An Example of Calculated Fitness Values 



Chromosome 


Fitness Value 


35782146 


144 


4382 7 561 


63 


6834 5721 


151 


74236815 


72 


21548637 


185 



20 



An example of selecting a pair of parent chromosomes from the population of Table 3 
based on their fitness might be to select 4 3 8 2 7 5 6 1, with a fitness of 63 and 7 4 2 3 6 
8 1 5 with a fitness of 72, supposing that 63 and 72 are the lowest fitness values of all the 
chromosomes in the population (even those not shown) and that a low fitness value is 
10 better. Alternatively, the fitness values could be maximized— it depends on what the 
fitness measures and the objective of the method. An example objective is to find a route 
through the locations that travels the minimum distance in the shortest time. An example 
of crossing over the pair, 4 3 8 2 7 5 6 1 and 7 4 2 3 6 8 1 5, at a randomly chosen point, 
the third position, to form offspring is shown in Table 4. To maintain a legal tour, i.e. an 
1 5 ordered sequence having no repeated numbers, the crossover operation may not 

reproduce a location in a child that is already in the child. For example, crossing over 4 3 
8 2 7 561 with 7 4 2 3 6 8 15 at the third position would create 4 3 8 3 6 8 1 5, if you 
started with 4 3 8 from the first parent and ended with 3 6 8 1 5 from the second parent. 
But 4 3 8 3 6 8 1 5 has duplicates of 3 and 8. In Table 4, rather than inserting duplicates, 
random numbers between 1-8 that do not already appear in the child are inserted in their 
place. In some embodiments, the starting point and ending point are unaltered in each 
iteration, e.g. the starting point may be close to the current location of transportation and 
the ending point may be close to a sorting location. These duplicates are shown as boxes 
and the inserted numbers are shown underlined in Table 4. 
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su 



10 



43827 561 



Table 4. An Example of Crossing over 

1 7 423681 5 



4 38D6D15 
7 42 OD5 6 1 



43826715 
43876215 
74 2 38561 
74 283 561 



An example of mutating offspring at each position on their chromosomes, with mutation 
probability p m between about 0.001 and 0.01 is no mutation at all. Alternatively, one or 
more chromosomes could be mutated at one or more positions while maintaining a legal 
tour. For example, 4 3 82671 5 could be mutated to483 2 6715, where 3 and 8 
have been swapped in the second and third positions. An example of a series of 
generations is shown in Table 5. In this example, the fitness values converged to paths 
starting with 4 3 8 and having fitness values between 34 and 40 by generation 3, with the 
best route being 4 3 8 2 6 7 1 5, with a fitness of 34. If the method stopped with 
generation 3, since the best fitness, 34, was found in both generation 2 and generation 3, 
then the solution route would be 4 3 82671 5. 



Table 5. An Example A Series of Generations 



Generation 1 


Fitness Values 


35782146 
43827 561 
6834 57 21 
74236815 
2154863 7 


144 

63 
151 
72 
185 


Generation 2 


Fitness Values 


438267 1 5 
4387621 5 
74238561 
74283561 


34 
65 
71 
73 


Generation 3 


Fitness Values 


438267 1 5 
43827615 
43827 651 


34 
32 
40 
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FIGS. 13 A and 13B show a more detailed flow chart of alternate embodiments 
1300 of the method 1200 of searching among locations for paths of FIG. 12. In one 
embodiment, each population of chromosomes has 2 k members, where the set of 
locations is of size k. In another embodiment, crossing over the pairs at the randomly 
chosen point 1208 is carried out by implementing a greedy crossover operation 1304. In 
another embodiment, the greedy crossover operation comprises generating offspring in 
the new population by 1306: placing all adjacencies shared by a first and a second parent 
in offspring 1308, alternating other adjacencies to be placed in offspring between the first 
and second parents 1310, and, when an adjacency produces a conflict, placing a random 
location in the offspring instead so that a legal tour is maintained 1312. In an alternate 
embodiment, the greedy crossover operation comprises generating offspring by 1314: 
repeatedly placing a shortest parental adjacency to a location not yet visited in an 
offspring, if such shortest parental adjacency exists 1316, otherwise placing a random 
adjacency in the offspring that maintains a legal tour 1318. An adjacency is the 
relationship of two locations next to each other in a chromosome, and in the real world, 
an adjacency represents a path, such as a road, between the locations, such as cities. For 
example, in the path a, b, c, h,f, d, e, g, locations c and h are adjacent. Also, h and/are 
adjacent. 

FIG. 14 shows another more detailed flow chart of alternate embodiments of the 
method of searching among locations for paths of FIG. 12. In one embodiment, selecting 
the pairs of parent chromosome from the current population based on their fitness 1206 is 
carried out by implementing a CHC selection operation 1402. In another embodiment, 
the CHC selection operation comprises 1404: placing the pairs of parent chromosomes in 
the new population 1406, merging the new population with the current population based 
on fitness 1408, and, upon determining premature convergence, retaining a number of the 
current population based on fitness and replacing the rest of the current population with 
randomly generated chromosomes 1410. Premature convergence occurs when relatively 
few generations are created before the fitness doesn't get much better in the next 
generation and the paths are very similar or identical with similar or identical fitness 
values. 
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FIG. 15 is another example embodiment of a method for solving a TSP. FIG. 15 
shows a flow chart of an example embodiment of a method 1500 for finding distribution 
chains, which are one part of a supply chain to satisfy customer demand to be 
programmed on a processor in a computer system. One aspect of the present invention is 

5 a computer system for finding distribution chains to satisfy customer demand, 
comprising: a storage device, an output device, and a processor programmed for 
repeatedly performing a method 1500. The computer system may be similar to the one 
shown in FIG. 19. The method 1500 comprises: determining current orders 1502, 
determining a set of locations to visit which comprises locations of suppliers with 

10 products available for the current orders 1504, determining a starting location in the set of 
locations 1506, determining a distribution center proximate to customer locations, the 
customer locations being delivery locations for each customer who placed one of the 
current orders 1508, applying a genetic algorithm to find a first route starting at the 
starting location, ending at the distribution center, and visiting each location in the set of 

15 locations, while minimizing distance and time 1510, providing the first route on the 
output device 1512, and storing any undelivered products in a buffer on the storage 
device to be included in the current orders and removing the orders that the first route 
will fill from the current orders 1514. 

FIG. 16 shows a more detailed flow chart of alternate embodiments 1600 of the 

20 method 1500 for finding distribution chains to satisfy customer demand of FIG. 15. In 
one embodiment, the method 1500 further comprises computing the distance between all 
locations, including each location in the set of locations, the starting location, and the 
distribution center 1602. In another embodiment, the method 1500 further comprises 
computing the distance between each customer location and the distribution center 1604. 

25 In another embodiment, the method 1500 further comprises applying the genetic 

algorithm to find a second route from the distribution center to each customer location 
that minimizes distance and time 1606. In another embodiment, the method 1500 further 
comprises providing the second route on the output device 1608. 

FIG. 17 shows example embodiments of data signals 2702 and 1708 to be 

30 communicated by an example embodiment of a method for solving the TSP. One aspect 
of the present invention is a data signal 1702 having a data structure stored thereon, 
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comprising: a first field containing locations 1704 and a second field containing a request 
for a route through the locations in the first field 1706. Another aspect of the present 
invention is a data signal 1708 having a data structure stored thereon, comprising: a field 
containing a route which comprises ordered locations to be visited. In one embodiment, 
the route of data signal 1708 is a delivery route. In another embodiment, the route of data 
signal 1708 is a pickup route. 

FIG. 18 shows an example embodiment of a method 1800 of communication for a 
user interface to a method for solving the TSP. One aspect of the present invention is a 
method 1800 of communication, comprising: receiving a plurality of locations and a 
request for a route through the locations 1802 and sending, in response to the request, the 
route 1804. This method may be implemented on a server computer in a client-server 
system. 

FIG. 19 shows an example of a suitable computing system environment 1900 for 
implementing embodiments of the present invention, such as those shown in FIGS. 6 A, 
6B, and 7-18. Various aspects of the present invention are implemented in software, 
which may be run in the environment shown in FIG. 19 or any other suitable computing 
environment. The present invention is operable in a number of other general purpose or 
special purpose computing environments. Some computing environments are personal 
computers, server computers, hand held devices, laptop devices, multiprocessors, 
microprocessors, set top boxes, programmable consumer electronics, network PCS, 
minicomputers, mainframe computers, distributed computing environments, and the like. 
The present invention may be implemented in part or in whole as computer-executable 
instructions, such as program modules that are executed by a computer. Generally, 
program modules include routines, programs, objects, components, data structures and 
the like to perform particular tasks or implement particular abstract data types. In a 
distributed computing environment, program modules may be located in local or remote 
storage devices. 

FIG. 19 shows a general computing device in the form of a computer 1910, which 
may include a processing unit 1902, memory 1904, removable storage 1912, and non- 
removable storage 1914. The memory 1904 may include volatile 1906 and non- volatile 
1908 memory. Computer 1910 may include — or have access to a computing 
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environment that includes — a variety of computer-readable media, such as volatile 1906 
and non- volatile memory 1908, removable 1912 and non-removable storage 1914. 
Computer storage includes RAM, ROM, EPROM & EEPROM, flash memory or other 
memory technologies, CD-ROM, digital versatile disks (DVD) or other optical disk 
storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic 
storage devices, or any other medium capable of storing computer-readable instructions. 
Computer 1910 may include — or have access to a computing environment that 
includes — input 1916, output 1918, and a communication connection 1920. The 
computer 1910 may operate in a networked environment using a communication 
connection 1920 to connect to one or more remote computers. The remote computer may 
include a personal computer, server, router, network PC, a peer device or other common 
network node, or the like. The communication connection 1920 may include a local area 
network (LAN), a wide area network (WAN) or other networks. Data signals, such as the 
ones shown in FIG. 18, may be communicated using the communication connection 
1920. A data signal may be sent over a wired network or direct- wired connection or 
wireless media such as acoustic, RF, infrared or other wireless media. 

It is to be understood that the above description is intended to be illustrative, and 
not restrictive. Many other embodiments will be apparent to those of skill in the art upon 
reviewing the above description. For instance, TBD. Further embodiment will be 
apparent to those of ordinary skill in the art. The scope of the invention should, 
therefore, be determined with reference to the appended claims, along with the full scope 
of equivalents to which such claims are entitled. 
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